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Abstract 

Model legumes such as Lotus japonicus have contributed significantly to the understanding of symbiotic 
nitrogen fixation. This insight is mainly a result of forward genetic screens followed by map-based cloning 
to identify causal alleles. The L. japonicus ecotype 'Gifu' was used as a common parent for inter-accession 
crosses to produce F2 mapping populations either with other L. japonicus ecotypes, MG-20 and Funakura, 
or with the related species L. fiiicaulis. These populations have all been used for genetic studies but seg- 
regation distortion, suppression of recombination, low polymorphism levels, and poor viability have also 
been observed. More recently, the diploid species L. burttii has been identified as a fertile crossing partner 
of L. japonicus. To assess its qualities in genetic linkage analysis and to enable quantitative trait locus 
(QTL) mapping for a wider range of traits in Lotus species, we have generated and genotyped a set of 
163 Gifu x L. burttii recombinant inbred lines (RILs). By direct comparisons of RIL and F2 population 
data, we show that L. burttii is a valid alternative to MG-20 as a Gifu mapping partner. In addition, we dem- 
onstrate the utility of the Gifu x L. burttii RILs in QTL mapping by identifying an N/r7-linked QTL for 
Sinorhizobium fredii nodulation. 
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1 . Introduction 

Symbiotic associations between plants and microor- 
ganisms have a large impact on ecosystems and 
plant production worldwide, and two prominent 
symbionts are phosphorus delivering fungi of the 
Glomeromycota, 1 forming arbuscular mycorrhiza, 
and nitrogen fixing soil bacteria, rhizobia. Since the 
most widely used model plant Arabidopsis tbaliana is 



unable to host these symbionts, two model legumes, 
Lotus japonicus and Medicago truncatula, were intro- 
duced to allow the genetic dissection of symbiotic 
processes. 

These model legumes are diploid, inbreeding species 
with relatively small sequenced genomes, 2,3 and they 
have allowed the isolation of a number of genes 
needed for mycorrhizal and rhizobial symbiosis. 4 
Most of these genes were isolated through map- 
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based cloning, but transposon and retrotransposon 
tagging have also been used. 5,6 

Selection of mapping parents is important in order 
to establish F2 populations that allow successful 
causal gene identification by map-based cloning. 
Ideally, the number of inter-parental polymorphisms 
should be high, while the F2 populations should be 
characterized by low levels of segregation distortion 
and suppression of recombination. For Lotus, the 
first genetic maps were produced using amplified 
fragment length polymorphism and microsatellite 
markers to genotype F2 offspring. The mapping popu- 
lations were established using L. japonicus Gifu as the 
common parent either in crosses with MG-20 and 
Funakura or in an interspecific cross with L. filicau- 
lis. 7 ~ 9 The level of polymorphism between 
L. filicaulis and Gifu is high but the F2 populations suf- 
fered from low viability and from severe segregation 
distortion in large genomic regions. 9,10 Whereas Gifu 
and Funakura were genetically very similar, a relative- 
ly high level of polymorphism was found between 
Gifu and MG-20, and the F2 populations showed 
good viability combined with little segregation distor- 
tion. 7 A dense linkage map with several hundred 
microsatellite markers was therefore developed, and 
most map-based cloning projects in Lotus have been 
based on crosses between Gifu and MG-20. 2 The 
major disadvantage of these populations is the large 
translocation between the top of Gifu chromosome 
1 and the bottom of MG-20 chromosome 2, 7 which 
corresponds to a genomic region with complete 
suppression of recombination. In parallel with these 
early genetic studies, Gifu-based recombinant 
inbred lines (RILs) were developed from crosses with 
MG-20, Funakura and L. filicaulis. 8,: 0,1 1 

To allow map-based cloning in the intervals display- 
ing the suppression of recombination in Gifu x MG-20 
populations, L. burttii originating from West Pakistan 
was introduced as a fourth crossing partner. 1 2,1 3 
Compared with MG-20 and L. filicaulis, L. burttii 
showed an intermediate level of polymorphism with 
respect to Gifu. This high marker density combined 
with good viability in offspring from crosses with Gifu 
and with the lack of suppression of recombination at 
the top of chromosome 1 made L. burttii a promising 
source of genetic material. 1 3 Despite these favourable 
characteristics, Gifu x L. burttii mapping populations 
have been sparsely used, probably because they have 
not been extensively characterized. This issue can be 
addressed by generating and analysing RILs, since they 
capture and stabilize a large number of recombination 
events, which can be exploited for the characterization 
of genome structure and recombination landscapes, as 
well as for quantitative trait locus (QTL) mapping. 

Here, we describe the development and genotyping 
of RILs derived from a Gifu x L. burttii F2 population. 



By analysing RIL data, we show that L. burttii is a 
valid alternative to MG-20 as a crossing partner to 
Gifu in mapping schemes, and we demonstrate that 
the Gifu x L. burttii RILs can be used to guide the se- 
lection of mapping parents. To illustrate the use of the 
RILs in QTL mapping, we investigated the nodulation 
of L. burttii by Sinorbizobium fredii HH1 03, which is a 
fast-growing rhizobial strain that nodulates soybean 
and a number of other legumes. 14,1 5 We found that 
L. burttii differs from Gifu and MG-20 in forming nitro- 
gen-fixing nodules with S. fredii HH103 and demon- 
strated the utility of the RILs by identifying a major 
Nfrl -linked QTL controlling S. fredii compatibility. 

2. Materials and methods 

2. 7. RIL propagation and genotyping 

Inbred lines of L. japonicus Gifu B-129 16,17 and 
L. burttii B-303 13 were used for crossing. Gifu x 
L. burttii RILs were developed from an F2 population 
by selfing until F8 and F9 using single-seed descent. 
A total of 1 63 lines have reached F8. 

Microsatellite markers developed to detect poly- 
morphism between Gifu and MG-20 2,18 were tested 
for polymorphism between Gifu and L. burttii (http:// 
www.kazusa.or.jp/lotus/markerdb_indexhtml). Additional 
markers were developed for sequenced transformation 
competent artificial chromosome and bacterial artifi- 
cial chromosome clones by searching for microsatellite 
sequences with the helpof WebSat. 19 Newly developed 
markers are listed in Supplementary Table S7. The 
markers were either analysed on 2-3% agarose gels 
or using an ABI3730 48 capillary sequencer (Applied 
Biosystems) in conjunction with the proprietary 
GeneMappersoftware (version 3.7) usingfluorescently 
labelled primers. One hundred and forty-six lines were 
scored with 9 7 microsatellite markers distributed 
across all six chromosomes (Supplementary Table S1 ). 

2.2. Analysis of RIL data and QTL mapping 

The RILgenotypingdata presented in Supplementary 
Tables S1 and S2 were imported into R/qtl 20 version 
1 .2 1 -2 (read.cross), converted to RIL format (conver- 
t2riself), and recombination fraction and logarithm 
of odds (LOD) scores plotted (plot.rf; Fig. 1A and B). 
For the QTL analysis, the genetic map was estimated 
(est.map and replace.map) and missing genotype 
values extrapolated (mqmaugment). QTL scans were 
then conducted using 1 000 permutations to deter- 
mine significance thresholds (mqmpermutation) and 
plotted (mqmplot.permutations). 

To allow the comparison of physical distances and 
recombination frequencies, genotyping markers 
were anchored to the L. japonicus MG-2 0 genome 
version 2.5 using BLAST. 21 Recombination fractions 
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Figure 1. Marker linkage and recombination landscapes for Gifu x MG-20 and Gifu x L. burttii RIL populations. (A and B) Recombination 
fractions (rf) and LOD scores for all marker pairs. Recombination fractions are in the upper left triangle and the LOD scores are in the 
lower right triangle. The LOD scores are for a test of r = 0.5. Red corresponds to a large LOD or a small recombination fraction, indicating 
a strong marker linkage, while blue is the reverse. Missing values appear in light grey. All plotted values are listed in Supplementary 
Tables S3 and S4. (A) Gifu x MG-20 RILs. (B) Gifu x L. burttii RILs. (C) The physical distance (in kb) corresponding to a 
recombination fraction of 1 % is plotted at each marker position. A high value indicates a low rate of recombination. (D) The fraction 
of L. japonicus Gifu genotypes at each marker position is plotted. A fraction of 0.5 (black line) is expected in the absence of 
segregation distortion. 
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were calculated using MapMaker3 22 with error detec- 
tion on. Physical distances per recombination fraction 
were then calculated (Supplementary Tables S5 and 
S6). To avoid inflated values during plotting 
(Fig. 1 C), a minimum recombination fraction of 
0.004 was used. 

2.3. Nodulation assays 

Seeds were germinated and plants grown on 1 /4 
B&D medium as described. 23 Sinorhizobium fredii 
HH1 03 was grown in YMB medium for 2 days and 
then diluted to an OD600 of 0.01-0.02 before in- 
oculating with 400 jutl per plant by pipetting the sus- 
pension directly onto the root. The nodulation 
phenotype was scored 6 weeks post-inoculation. 

3. Results and discussion 

3. 1 . Generation and genotyping of 1 63 Gifu x L. burttii 
RILs 

Characterization of populations derived from crosses 
between Gifu and L. burtti is required to promote the 
exploitation of the high level of Gifu//.. burttii poly- 
morphism for genetic mapping. For this purpose, and 
to provide the legume community with a stable 
genetic resource for QTL and recombination rate ana- 
lysis, Gifu x L. burttii RILs were generated by selfing. 
Starting from 200 F2 plants, 163 lines have reached 
generations 8 or 9. One hundred and forty-six of 
these lines were genotyped by scoring 97 microsatel- 
lite markers (Supplementary Table S1), 2 and a Gifu x 
L burttii genetic map, syntenic to the Gifu x MG-20 
map, 2,7 was constructed (Supplementary Table S1). 
The RILs have been deposited at LegumeBase, 
Miyazaki, Japan, where they are available for the scien- 
tific community (http://www.legumebase.brc.miyazaki- 
u.ac.jp/). The genotypes of the Gifu x MG-20 RILs can 
be found at ftp://ftp.kazusa.or.jp/pub/lotus/Rlline/ 
and in Supplementary Table S2. 

3.2. Gifu x L. burttii and Gifu x MG-20 mapping 
crosses yield complementary information 

To assess whether Gifu x L. burttii populations 
would offer genetic resolution in chromosomal 
regions where the suppression of recombination is 
prominent in populations derived from crosses 
between Gifu and MG-20, we compared recombin- 
ation landscapes based on RIL genotyping data. 
Analysing the pairwise recombination fractions and 
LOD scores for all marker combinations in the Gifu 
x MG-20 RILs resulted in a clear marker linkage 
signal at the previously described Gifu/MG-20 trans- 
location between chromosomes 1 and 2 (Fig. 1A 
and Supplementary Table S3). 7 This signal was 
absent in the Gifu x L. burttii RIL data (Fig. 1 B and 



Supplementary Table S4), confirming previous obser- 
vations. 13 However, markers at the top of chromo- 
some 1 and the bottom of chromosome 4 appeared 
linked in the Gifu x L burttii population, suggesting 
an altered L. burttii chromosomal structure in this 
region (Fig. 1 B and Supplementary Table S4). 

To investigate if these differences in the genome 
structure were reflected in differential recombination 
rates, the genotyping markers were anchored on the 
MG-20 release 2.5 genome assembly, and physical 
distance per recombination fraction was calculated 
and plotted (Fig. 1C and Supplementary Tables S5 
and S6). 

Both a relatively high inter-marker distance and an 
incomplete genome assembly contribute to the un- 
certainty of this analysis, and conclusions should be 
drawn with care. We did, however, observe a higher 
recombination rate in the Gifu x L. burttii population 
than in the Gifu x MG20 population (Fig. 1C) in the 
genomic interval located at the top of chromosome 
1 in the MG-2 0 reference sequence, as noted previ- 
ously. 13 A higher recombination frequency for the 
Gifu x L burttii population was also found in the 
symW region near the top of chromosome 4. On 
the other hand, this population appeared to suffer 
from suppression of recombination in the central 
part of chromosome 1, but not in the chromosome 
1 and 4 regions displaying the marker linkage. 
Overall, it appears that the two populations show 
comparable levels of recombination and that the re- 
combination frequency varies significantly with both 
genomic location and population. In addition, the 
Gifu x L burttii cross showed a balanced segregation 
of genotypes, except for the central region of chromo- 
some 2 (Fig. 1 D). In conclusion, L. burttii appears to 
be a valid alternative to MG-20 as a Gifu crossing 
partner for generating Lotus mapping populations. 

3.3. Mapping of tbe symbiotic mutant sym 1 0 

demonstrates the utility of RIL-guided selection 
of mapping parents 

Fine-mapping of the symlO locus located on 
chromosome 4 near 1 0 cM 10 has proved challenging 
using a Gifu x MG-20 population, as no recombin- 
ation events could be detected between the flanking 
markers TM1 237 and TM0326 in a total of 977 indi- 
viduals. This limited the mapping resolution to 
~600 kb. The Gifu x L. burttii RILs documented here 
represent a genetic resource, which can be used to 
test the applicability of a mapping cross in a given 
genomic region. For the sym 10 region, the high local 
rate of recombination indicated by the RIL data sug- 
gested that a Gifu x L. burttii population could be 
used to increase genetic mapping resolution (Fig. 1 C). 

Motivated by this result, we generated and geno- 
typed a symWxL. burttii mapping population. 
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Figure 2. sym 7 0 mapping resolution with crosses to L. japonicus MG-20 and L. burttii. The number of informative recombinants is shown in 
the cross(es) where the marker is polymorphic. The delimited size of the syml 0 region is shown for the two crosses individually and with 
the combined mapping information. 



Although 977 individuals from the symW x MG-20 
population allowed us to narrow down the candidate 
region to ~600 kb, the 696 individuals from the 
sym 10 x L. burttii population provided additional in- 
formation, limiting the candidate region to ~300 kb 
(Fig. 2). 

In this case, the high recombination rate in the 
region of interest suggested by the RIL data was 
reflected in the F2 mapping population, leading to 
increased genetic resolution when combining infor- 
mation from the two populations. We propose that 
taking advantage of the RIL resources to probe recom- 
bination frequencies in genomic regions of interest 
prior to the selection of the mapping partner for F2 
fine-mapping populations could allow faster 
mapping progress in a number of cases. 

High-throughput sequencing of bulked mapping 
populations generated either by interspecies or by 
back-crosses has proved effective in identifying causal 
ethyl methanesulfonate (EMS) mutations. 24,25 For 
the sequencing analysis of interspecies crosses, 
marker density and recombination frequencies 
remain important parameters, and selecting the 
optimal mapping population cross is still important. 
Deep sequencing back-crossing methods also require 
genetic mapping information, and the induced EMS 
mutations have been used as markers in place of inter- 
species polymorphisms. 25 A number of interesting 
symbiotic mutants in Lotus are not generated by EMS 
mutagenesis but rather by unknown genomic 
changes induced by tissue culture regeneration, 26 
and accurate genetic mapping information remains 
critical for identifying these causal mutations. 

3.4. A major QTL near Nfr1 controls S. fredii 
compatibility 

In addition to facilitating the choice of mapping 
crosses for the identification of causal mutations, 




Figure 3. S. fredii infection phenotypes for L. burttii and Gifu. 
Six-week-old plate-grown plants inoculated with S. fredii. (A) 
L. burttii. (B) L. japonicus Gifu. Arrow points to red root nodule. 



the genotyped Gifu x L. burttii RILs should also be 
well suited for exploiting Lotus natural variation to 
identify genetic determinants of a wide range of 
traits through QTL analysis. 

Since we had observed a clear difference in nodula- 
tion with S. fredii HH1 03 between L. burttii and Gifu 
(Fig. 3), we chose this trait as an example for QTL ana- 
lysis and proceeded to record the S. fredii nodulation 
phenotype of all RILs. Twenty plants were scored per 
line, and the trait was quantified as the 'average 
number of pink nodules per plant'. 

Analysing the data using R/qtl initially revealed two 
significant peaks. One was present near the marker 
TM0002, which is closely linked to the Nod factor re- 
ceptor Nfrl , 27 and another was found at the bottom 
of chromosome 4 near marker TM0617 (Fig. 4A). 
Further, a minor peak was seen at the top of chromo- 
some 5. Rerunning the QTL analysis including the 
markers corresponding to the three highest peaks as 
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Figure 4. QTL analysis of S. /red/7 nodulation. The trait analysed is 
the average number of nodules on Gifu x L. burttii RILs 
inoculated with S. fredii. (A) QTL analysis without cofactors. (B) 
QTL analysis including markers TM0002, TM061 7, and 
TM0494, corresponding to the three significant peaks in (A), as 
cofactors. Five and 1 0% false discovery rate levels are indicated 
by the upper and lower dashed lines, respectively. Light grey 
area delimits the results of 1 000 QTL analyses on permuted 
data, which were used to determine the significance thresholds. 
Dashed vertical grey lines indicate chromosome ends. 

cofactors eliminated the peak on chromosome 4 
(Fig. 4B), indicating that it was an artefact caused by 
marker linkage between the top of chromosome 1 
and the bottom of chromosome 4 (Fig. 1 B). The vast 
majority of the phenotypic variation, therefore, 
appeared to be explained by the QTL on chromosome 
1 near Nfr1 . 

In conclusion, the Gifu x L. burttii RILs have allowed 
us to identify an Nfrl -linked locus, which supports an 



unusual rhizobial-/_otus symbiotic interaction. In con- 
trast to the N/r5-dependent host-range extension 
shown for L. filicaulis previously, 28 the Nfr5 locus 
does not seem to influence the interaction between 
L burttii and S. fredii. 

The present set of Gifu x L burttii RILs expands 
both the genotypic and the phenotypic space avail- 
able to Lotus researchers and promises to deliver 
new insight into the genetic mechanisms governing 
natural variation in rhizobial-legume interactions. In 
addition, the RIL population displays variation in 
flower colour, leaf shape, flowering time, and pod 
shattering, making it well suited for QTL mapping of 
these traits. 
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